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The universe can be made flat and smooth by undergoing a phase of ultra-slow (ekpyrotic) contrac- 
tion, a condition achievable with a single, canonical scalar field and conventional general relativity. 
It has been argued, though, that generating scale-invariant density perturbations, requires at least 
two scalar fields and a two-step process that first produces entropy fluctuations and then converts 
them to curvature perturbations. In this paper, we identify a loophole in the argument and intro- 
duce an ekpyrotic model based on a single, canonical scalar field that generates nearly scale-invariant 
curvature fluctuations through a purely "adiabatic mechanism" in which the background evolution 
is a dynamical attractor. The resulting spectrum can be slightly red with distinctive non-gaussian 
fluctuations. 



Two mechanisms are known for making the universe 
flat and homogeneous in accord with cosmological obser- 
vations. The first is inflation [T], a period of accelerated 
expansion occurring after the big bang during which the 
equation of state w is less than —1/3. The second is 
ekpyrosis [H [3] , a period of ultra-slow contraction with 
w I occurring before a bounce to an expanding phase. 
Both equations of state can be obtained with a single, 
canonical scalar field evolving along a potential and con- 
ventional general relativity. To obtain a nearly constant 
w, the potential must be positive and flat for inflation 
and steep and negative for ekpyrosis. In both cases, the 
scalar field and Newtonian gravitational potential obtain 
nearly scale-invariant fluctuations. 

In addition, the scalar field fluctuations automatically 
induce nearly scale-invariant curvature (density) pertur- 
bations in inflationary models. The same is not true for 
ekpyrosis. It has been argued that the spectrum for the 
curvature perturbation £ in ekpyrotic models with a sin- 
gle canonical scalar field is necessarily blue [3], incon- 
sistent with observations of the cosmic microwave back- 
ground and large-scale structure. A scale-invariant spec- 
trum can be obtained via an "entropic mechanism" [5], 
but this requires at least two scalar fields and a two-step 
process that first produces entropy fluctuations and then 
converts them to £ [6] . 

Another feature of inflation is that the cosmological 
background is a dynamical attractor, as can be seen from 
the behavior of £ « 8a/ a in the long- wavelength (fc — > 
0) limit. Since £, which measures the perturbation in 
the expansion history between distant Hubble patches, 
approaches a constant in this limit, the perturbation can 
be absorbed through a spatial diffeomorphism: hence, 
the background solution a(t) is an attractor. 

When compared to the known alternative cosmologi- 
cal models, it seemed that inflation is unique in having 
the combined properties of scale-invariance and dynami- 
cal attractor. For instance, the background is an attrac- 
tor and £ is constant at long wavelengths in single-field 
ekpyrosis 013], but, as mentioned earlier, its spectrum 
is blue, rather than scale-invariant. As another case in 



point, it is well-known that the equation for £ in a con- 
tracting, matter-dominated universe is identical to that 
of inflation, and hence results in a scale-invariant spec- 
trum [7j; but the growing mode for £ increases in time 
at long wavelengths, indicating that the background is 
unstable. This instability also plagues contracting mech- 
anisms relying on a time-dependent sound speed [8J. 

In this paper, though, we present a counterexample: a 
model based on a single, canonical scalar field in which 
the contracting cosmological background is a dynamical 
attractor and a scale-invariant spectrum of curvature per- 
turbations is generated via a purely "adiabatic mecha- 
nism." More precisely, the power spectrum is slightly 
red, I — n s = 0(1%), with distinctive /c-dependent non- 
gaussianity (compared to the entropic mechanism or in- 
flation). As with the entropic mechanism, this new adia- 
batic mechanism is observationally viable assuming that 
the bounce from contraction to expansion preserves the 
long-wavelength spectrum for Q. 

The counterexample relies on relaxing the usual as- 
sumption that w is nearly constant. The effect can be 
illustrated for fairly simple scalar potentials, such as 

V{<t>) = V {l-e-^), (1) 

with c 3> 1. (Here and throughout, we use reduced 
Planck units for which the gravitational constant is G = 
I/87T.) The regime of interest is the transition when the 
equation of state rises rapidly from e = |(1 + w) « 0, 
where the constant term dominates, to e « c 2 /2, where 
the negative exponential term dominates. The transi- 
tion, which occurs when ~ c" 1 logc and V is positive, 
produces an ekpyrotic phase during which scale- invariant 
adiabatic curvature perturbations are generated. 

The form for V(4>) is similar to inflationary examples, 
but here we are assuming a scenario in which the universe 
is slowly contracting just prior to the onset of the adia- 
batic mechanism. The transition in the equation of state 
occurs over a tiny range Acf> <C 1 around = 0, so (JID 
need only be a good approximation over this range and 
there is considerable freedom in the overall form of the 
potential 0. 
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Note that we focus in this paper on the finite period 
during which the adiabatic mechanism generates pertur- 
bations presuming some pre-existing contracting phase 
and some subsequent bounce and reheating. The various 
ways in which adiabatic mechanism can be embedded in 
a larger scenario are left for future work. 

Until now, all ekpyrotic potentials required V < in 
order to obtain the desired properties for the equation of 
state: w > 1 and nearly constant. For example, for the 
negative exponential potential alone, V{(f>) = — Vbe _c ^, 
the solution has H = 2/c 2 t and H = —2/c 2 t 2 , where 
time t — > 0~ corresponds to the bounce. Then, 



2 V 



H 
IP 



(2) 



a constant value much greater than unity for c » 1. 

By contrast, the scale-invariant curvature perturba- 
tions described in this paper are generated during the 
transition from the positive constant term dominating to 
the exponential dominating, which we term the "transi- 
tion phase" (to distinguish from the conventional ekpy- 
rotic phase that occurs subsequently). During this tran- 
sition, we will find that the equation of state parameter e 
grows as 1 jt 2 , obtaining the large positive values needed 
to smooth and flatten the universe. This behavior for e 
leads to a range of wavelengths for which the spectrum 
of curvature perturbations is scale invariant. 

As in inflation, the adiabatic mechanism is based on 
the equation of motion for the Fourier modes £fc with 
comoving wavenumbers fc; with a change of variable 
Cai( t ) = v k/z, the equation of motion takes the form 



k K = , 

z 



(3) 



where primes denote derivatives with respect to confor- 
mal time r, and 



z = a(T)V2e(r) 



(4) 



In the "transition phase" , the scale factor is nearly con- 
stant (a(r) w 1) — this is the slowly contracting back- 
ground typical of ekpyrotic cosmology. It follows that 
conformal time and cosmological time are approximately 
the same: t « r. Meanwhile, as we will see, the equation 
of state parameter satisfies e ~ 1/t 2 s» 1/r 2 . Therefore, 
^tran(T) ~ ( — T ) — exactly as in inflation, where e is 
nearly constant and a(r) w r)! The spectrum is 
therefore scale-invariant. Note, morever, that the grow- 
ing mode solution for k — ¥ is ( — v/z — > constant: the 
"transition phase" evolution is a dynamical attractor. 

The growing mode solution is oscillatory for 
ky/\z/z"\ > 1 and constant for ky/\z/z"\ -C 1 if z" / z > 
0; the cross-over from oscillatory to constant behavior, 
called freeze-out, occurs when the wavelength of the mode 
is of order H^ czc = y/\z/z"\, the freeze-out horizon ra- 
dius. The notation is chosen to remind the reader that 
the freeze-out horizon, the relevant scale for analyzing 



perturbations, is of order the Hubble horizon H^ 1 for 
constant e; in this paper, though, it will be important to 
distinguish the two horizons. 

The conditions e ~ 1/t 2 arises naturally as a back- 
ground solution for simple potentials, such as 0. As 
usual in ekpyrotic cosmology, the background scalar field 
evolution is insensitive to the slowly-contracting back- 
ground and is therefore driven by V ((f)): 



4> ~ — cVge 



(5) 



We will check a posteriori that gravity can indeed be 
neglected in this equation. Note that the evolution of <f> 
is oblivious to the constant term in V(4>). The solution 
is therefore of the standard ekpyrotic form: 
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Substituting into H — —(p 2 /2 gives H s» —2/c 2 t 2 , which 
can be immediately integrated: H(t) = 2/c 2 t + Hq. At 
sufficiently early times, H is nearly constant, with the 
constant H fixed by the Friedmann equation: 3Hq w 
4> 2 /2 + V(<f>) — Vq. This defines the transition phase. 
In other words, the solution ^ is such that the kinetic 
energy nearly cancels the exponential term in , leaving 
the constant Vq term as the dominant contribution to H: 



-Htran(i) 
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The transition phase ends when the constant term no 
longer dominates, which can be read off from ([7]): 



'•end— tran 
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For t > fend- train we have H(t) ss 2/c 2 t, and the solution 
reduces to a standard ekpyrotic scaling phase. 
Meanwhile, the equation of state follows from 
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Deep in the transition phase, \t\ S> |icnd-tran|) wc have 
e ~ 1/t 2 , as desired to generate a scale- invariant spec- 
trum for £. For \t\ <C |t e nd-tran|j however, we have 
e w c 2 /2, indicating again an ekpyrotic scaling phase. 

Let us check the validity of assuming H(f> <C cV^e^ c< ^ 
in (jil. Substituting H m —y/Vo/Z and the solution 
it is easily seen that this approximation is consistent for 



* > t h , 



eg— tran 



(10) 



The transition is shorter than a Hubble time, and hence 
the scale factor, a tra n(t) ~ 1 — y/Vo/3 t is approximately 
constant throughout — the universe is nearly static. 
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Although H is nearly constant during the transition 
phase, this is emphatically not a contracting de Sitter 
universe. There are many equivalent ways to see this. 
First, as was already mentioned, the transition from e rs 
to e « c 2 /2 occurs within a Hubble time, and hence 
a(t) ~ 1. Moreover, the rate of change of e is never 
small: T) = .H" _1 dlne/di ranges from 0(1) to c 2 3> 1. 

The transition phase is followed by a standard scaling 
ekpyrotic phase in which e w c 2 /2. For |t| <C |i e nd-tran|, 
the background is therefore slowly-contracting with [2] 



^scaling (^) ( ^) 
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This standard ekpyrotic phase is part of the story as well: 
it enables the Hubble horizon to shrink while preserv- 
ing the scale-invariant spectrum such that the modes are 
well-outside the Hubble horizon before the bounce. 
As for C, we have from ^ and a(t) ~ 1 that 



^■tran (<) =a(i)v / 2^t) 



2v/3 
/VoH) 



(12) 



within the transition phase when \t\ 3> |t e nd-tran|- Thus, 
z" / z rj z/z « 2/i 2 , and the solution to ^ for the mode 
functions, with the adiabatic vacuum choice, is 
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These mode functions are identical to those in de Sit- 
ter space, except that their amplification relies here on a 
rapidly-changing e(t), not a rapidly-changing a(t). The 
long-wavelength amplitude for the curvature perturba- 
tion, = Vk/z, is therefore scale-invariant: 



(14) 



The range of these modes is set by the duration of 
the transition phase: fc max /fc min = tbeg-tranAend-tran 
c 2 /2. 

The transition phase is an attractor. Following [4], 
we can compute perturbations to in synchronous 
gauge and show that they all go to zero as 0(fc|t|) or 
as log(t/f e nd-tran), except for the C term which goes to 
a constant as usual. Thus the metric approaches its un- 
perturbed form. Similarly, 5(f) = 0(k\t\). 

In terms of local observables, the 3-curvature scalar 
tends to a constant, ~ k 2 £/a 2 , as usual. Meanwhile, 
the extrinsic curvature tensor, K t j = gij/2, receives a 
growing contribution from the log term in g^-: 
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However, recall that the background itself evolves away 
from de Sitter as H(t) = -y/V /3 + 2/c 2 t. Thinking of 
this as a correction over K = —y/SVo, we find 
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FIG. 1: Numerical computation of the perturbation ampli- 
tude fc 3 / 2 ^ vs. k generated by the adiabatic mechanism. The 
behavior of modes with larger and smaller k depends on the 
larger scenario in which the mechanism is embedded and be- 
yond the consideration of this paper. 



That is, the perturbation (15) behaves identically to the 
backg rou nd. Even more convincing is the observation 
that (15) is less singular than SH/H ~ l/< 2 resulting 



from a constant time shift of the background solution. 

To incorporate our mechanism in an observationally 
acceptable model, there are three further considerations: 
Hubble horizon exit: Although modes are amplified dur- 
ing the transition phase, they have not yet exited the 
actual Hubble horizon. Since both a(t) and H{t) are 
constant throughout this phase, modes remain within the 
Hubble horizon. Mode amplification nevertheless occurs 
because the freeze-out horizon radius, Hz} = y/\z/z"\, 
shrinks, thanks to the rapidly-changing e(t). 

The exit of modes outside the actual Hubble horizon is 
achieved during the subsequent scaling ekpyrotic phase, 
since |-ff| -1 ~ \t\ shrinks in this period. But since a(t) ~ 
(— t) 2 / c with c> 1 during this phase, the background is 
very slowly contracting, and the scale-invariant spectrum 
is preserved in the process. In other words, instead of 
= Ck ~ fc _1 / 2 z _1 , as is the case when the wavelength is first 



well below H^ CCZQ , the appropriate initial conditions for 
the mode functions at the onset of the ekpyrotic phase is 
modified: (k ~ /c _3 / 2 z _1 . Matching at horizon crossing 
gives A: 3 / 2 ^. ~ fc -1 ^ 1 " 6 ), which is nearly scale-invariant 



for 



72 > 1. 



These conclusions are borne out by numerical analy- 
sis. Using z — 2c(—t) 2 / c /(l + c 2 H t) to cover the tran- 
sition and scaling phases, we integrate Q with c = 140 
and \H \ = 5 • 1(T 4 over the interval -2.5|iJ r 1 < t < 
-5 ■ Kr 10 |iJ o r\ for the range of modes 0.02|iJ | < k < 
2 • 10 4 |i?o|- The shortest-wavelength mode is therefore 
barely outside the Hubble radius by the end of the inte- 
gration, which occurs deep in the ekpyrotic scaling phase. 
Figure [T] shows the resulting spectrum, with range of 
scale-invariant modes spanning a factor of 10 4 ~ c 2 . 

We can immediately derive constraints on Vq and c. 
Since the scale-invariant modes are on scales smaller than 
« -73/Vb, a necessary condition is that the co- 
moving scale 1 / vVq encompasses the entire observable 
universe. If the value of H at the end of the ekpy- 
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rotic phase, -ff c k-cnd, is comparable to that at the onset 
of the expanding, radiation-dominated phase, we must 
have -ffck-cnd/v^ Z ^ek— end-^ek— end/V^today-^today) ^ 
yj -ffck-cnd / H to day , where we have assumed a radiation- 
dominated evolution until the present epoch for simplic- 
ity. In other words, 



< 



\J -ffck-end-ff' 



today 



30 rrl/2 



(17) 



For Grand-Unified (GUT) reheating scale, -ff e k-cnd 
10 12 GeV, and the above condition is satisfied for Vq 



TeV 4 . Setting C 



from ( 14 ) that c 



■> 10~ 5 on large scales then requires 
10 28 . For electroweak (EW) reheat- 
ing, _ff e k-cnd ~ meV, we similarly get Vq ~ MeV 4 and 
c = 10 40 . 

Although we have focused on pure exponential poten- 
tials, for simplicity, the exponentially large values of c re- 
quired can be achieved effectively in the Conlon-Quevedo 
potential [ID] . V{4>) ~ exp(— 70 4 / 3 ), for large 4>. 
Spectral Index: Since the values of c of interest are expo- 
nentially large, the departures from scale-invariance are 
thus far unobservable. However we can generate a small 
red tilt, as favored by observations [TT|, by generaliz- 
ing ([!]) to allow for a slowly- varying exponent c((f>): 



V(cj>) = V + U(^) =V (l- e~S** 



(18) 



In the transition-phase approximation, the Hubble pa- 
rameter is constant, H s» — y/Vo/3, whereas U((f>) cancels 
against the kinetic energy of (f>. Therefore, 
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Since a(t) ps 1 as before, it follows that z ~ y/—U((f)). Us- 
ing the evolution equation <f> = — \f— 2U (</>), it is straight- 
forward to show that z" / z ps z/z ps —U^, where 
can be shown to satisfy ps — (2/< 2 )(l + 3c ^/c 2 ) [5]. 
Substituting in ([3|, the resulting tilt is 

n s - 1 = -4^| . (20) 

Since is decreasing in our solution, the spectral tilt will 
be slightly red if c ^ > 0. For instance, if c{4>) changes 
smoothly by 0(c) during the transition, then n s — 1 ss 
-4(Ac)/(c 2 A</») ps -4/(cA0). Using |6|, we obtain 



(21) 



log(^b e g_tran/ ^end— tran ) log(c 2 /2) 



For 10 40 > c > 10 28 , ranging from EW to GUT-scale 
reheating, we obtain n s ps 0.98. Since we dropped 0(1) 
factors, the generic answer is 1 — n s ps few %, in good 
agreement with observations |llj . 

Non-Gaussianity: As with the entropic mechanism [12] . a 
steep non-linear potential generates large non-gaussianity 
compared to simple inflationary models. In this case, the 
steepness increases during the transition phase, so modes 



that freeze out later have larger non-gaussianity; i.e., the 
non-gaussianity is fc-dependent. The precise calculation 
is lengthy and will be presented in [5J. The result is 

(CfciCteCfcs) = 4(27r) 3 5(fc! + k 2 + fc3)fc 6 |C*|M/IL *? with 
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where fc = fci + &2 + &3- This amplitude peaks for 
equilateral configurations (fci w fc 2 ~ fe) and vanishes 
in the squeezed or local limit (£3 <C fci ps fc 2 ) [SJ. 
As usual, the magnitude of „4 can be characterized by 
30*Afc, = fc/ 3 /£; 3 , which in our case gives 
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Note that /^ U1 increases with k and is of order c 2 for 
the shortest-wavelength mode generated during the tran- 
sition phase (& max = 1/Wd-tran ~ c 2 /VVo)- Given 
that exponentially large values of c are needed to have 
£ ~ 10~ 5 , there is the danger that the non-gaussianity 
grows unacceptably large for some range of k. (Ref. [13] 
recently pointed out additional problems with the pure 
exponential model.) The problems arise because the 
ekpyrotic phase with large c is maintained longer than 
needed. The simple solution is to alter the pure expo- 
nential potential so as to terminate the transition phase 
before the modes with unacceptably large /nT"' are 6 en ~ 
erated; for example, if c(<j>) decreases smoothly to some 
b <C c after the transition phase has generated an ob- 
servationally acceptable range of scale-invariant fluctu- 
ations. Current observations require £ ~ 10~ 5 over a 
range of k spanning a factor of 10 3 , which can be achieved 
while keeping the non-gaussianity acceptably small. Ter- 
minating the transition phase results in the perturbation 
amplitude being suppressed on smaller scales: the spec- 
trum tilts strongly to the red and then flattens out at 
an exponentially smaller amplitude with an acceptable 
non-gaussianity (/nlC ^ 1) throughout. The pertur- 
bative regime is consistently valid, both classically and 
quantum mechanically, at all times in this model [3J. 

We close by noting that the generation of a scale- 
invariant spectrum of adiabatic curvature perturbations 
over a range of k spanning 10 3 or so was once a theoreti- 
cal speculation and, after the introduction of inflationary 
cosmology, a theoretical prediction. Today it is a rigorous 
criterion for a viable cosmological model, given the ob- 
servations of the microwave background and large-scale 
structure. What has been learned from the ekpyrotic 
model is that the spectrum is not unique to inflation. 
With this paper, we have introduced an adiabatic ekpy- 
rotic mechanism that generates scale-invariant adiabatic 
perturbations, within a background that is a dynamical 
attractor. This mechanism also predicts non-gaussianity 
and a spectrum of gravitational waves that is observa- 
tionally distinguishable from inflation [3J. 
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